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The feed-forward relationship between diagnosis and prognosis is the foundation of both 
aircraft structural health management and the digital twin concept. Measurements of 
structural response are obtained either in-situ with mounted sensor networks or offline 
using more traditional techniques (e.g., nondestructive evaluation). Diagnosis algorithms 
process this information to detect and quantify damage and then feed this data forward to 
a prognostic framework. A prognosis of the structure’s future operational readiness (e.g., 
remaining useful life or residual strength) is then made and is used to inform mission- 
critical decision-making. Years of research have been devoted to improving the elements 
of this process, but the process itself has not changed significantly. Here, a new approach 
is proposed in which prognosis information is not only fed forward for decision-making, 
but it is also fed back to the forthcoming diagnosis. In this way, diagnosis algorithms can 
take advantage of a priori information about the expected state of health, rather than 
operating in an uninformed condition. As a feasibility test, a diagnosis-prognosis feedback 
loop of this manner is demonstrated. The approach is applied to a numerical example 
in which fatigue crack growth is simulated in a simple aluminum alloy test specimen. A 
prognosis was derived from a set of diagnoses which provided feedback to a subsequent set 
of diagnoses. Improvements in accuracy and a reduction in uncertainty in the prognosis- 
informed diagnoses were observed when compared with an uninformed diagnostic approach. 


I. Introduction 


In the context of structural health management (SHM), digital twin focuses on the coupling of high- 
fidelity life predictions with in-situ damage diagnostics.’'? These two concepts are often considered as 
separate entities. Damage diagnosis is an active area of research in which changes in modal properties, 
strain response, or guided wave scattering patterns are used to detect and quantify flaws or defects in a 
structure. Damage diagnosis is a process that involves multiple scales (i-e., damage is manifested at the 
microstructural level, and measurements are often taken at the macroscale, sometimes at large distances 
away). Uncertainties are introduced to the problem at each scale and via models used to form the diagnoses; 
therefore, quantifying these uncertainties is paramount if the gathered data is to be used for decision-making.? 
Damage prognosis is typically a two-step process that can also encompass multiple scales and suffers from 
multi-source uncertainty. To form a prognosis, the uncertainties in both the set of available diagnoses and 
the parameters of a chosen predictive model must be quantified. These uncertainties are then propagated 
through the model to produce a probabilistic prediction of remaining useful life (RUL).° 

Significant work has been invested in developing the relationship between damage diagnostics and prog- 
nostics. For example, Sankararaman et al. developed a prognosis framework to interface with damage 
diagnosis data and used it to predict fatigue crack growth in an aluminum cylinder with quantified uncer- 
tainties.° Peng et al. used Markov chain Monte Carlo (MCMC) sampling methods to combine diagnosis, 
prognosis and uncertainty quantification for an aluminum lap joint.© Similar methods have been applied 
to composite materials as well, as demonstrated by the particle filtering approach used by Chiachio et al.” 
to predict fatigue delamination growth. These are only a small sampling from the field, but all of these 
works share a common thread in that the relationship between diagnosis and prognosis is unidirectional. 
Furthermore, most approaches tend to treat the diagnostic problem as deterministic. 
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As uncertainty quantification becomes more prevalent in damage diagnostics research, an opportunity 
arises in which the current linear relationship between diagnosis and prognosis can be evolved into a more 
complex, cyclical relationship. More precisely, a feedback loop can be developed with prognosis feeding 
information about expected future damage conditions back to the diagnosis framework, as shown in Figure 
1. If the uncertainty quantification method used for diagnosis is based on a Bayesian ideology, then this 
information would take the form of a prior distribution. In situations where noise levels or lack of data 
inhibit Bayesian inference of the damage parameters, an informed prior distribution can allow for successful 
damage diagnosis. It should be noted that the overarching feedback concept could be extended to a more 
informed, experimental diagnostic approach as well. For example, when using guided wave-based diagnostics, 
the transmission and reflectivity ratios are highly dependent on the type and size of damage as well as the 
frequency of the impending wave.®:? A prognosis of the damage state at the time the next measurement 
is to be taken could allow for tuning of the excitation frequencies and a potential improvement of image 
resolution or probability of detection. The work presented herein focuses primarily on the idea of stochastic 
feedback, and experimental applications such as this are left for future work. 


Traditional, unidirectional relationship between diagnosis and prognosis 


Operational 
Decision 


Present study 


Proposed feedback loop 


Figure 1. Diagram of the proposed feedback loop. While feedback can be applied to measurement techniques as well 
as the diagnostic framework, the present work is only focused on providing stochastic feedback to the diagnosis. 


The hypothesis for the present study was formulated as follows: prognosis of the damage state at a future 
time instant can inform a future diagnosis made at that same time instant, resulting in an improvement of 
accuracy and a reduction in uncertainty. The feedback was assumed to be stochastic in nature, taking 
the form of a prior distribution fed to the Bayesian diagnostic framework. The work was intended to be a 
feasibility study; as such, to control for excessive randomness from model or experimental errors, a numerical 
fatigue crack growth simulation was conducted using finite element analysis (FEA). A set of diagnoses (i.e., 
strain-based estimates of crack length over time) were used to produce a probabilistic estimate of the crack 
length at future time instants corresponding to planned diagnoses. These predicted distributions of crack 
lengths were fed to the diagnostic framework as prior distributions, and the final diagnoses were made. To 
test the effect of the feedback loop, a diagnosis was also made with a noninformative prior distribution at 
each examined time. The accuracy and uncertainty of the two results were then compared relative to the 
ground truth simulation. 


II. Formulation 


Let the useful life of a structure be discretized by the times t;, where i = 1,...,p,...AOL. Damage 
prognosis involves making a diagnosis at each index up to tp, the time at which a prediction will be made. At 
this time, the diagnosis information for all past time steps are gathered and used to quantify the uncertainties 
in both the measurements and the predictive model. These uncertainties are propagated forward in time 
to some end of life (EOL), or the time at which failure is deemed to have occurred, tgor,. This concept is 
integral to both digital twin and SHM, as discussed in the introduction. However, predictive information 
from the time between t, and tgoz is often unused or discarded. Here, it is proposed that this information 
be utilized to inform the next damage diagnosis measurement at time t, + 7, where T is an arbitrary step 
forward in time. For any 7, a distribution of predicted damage states is available that can be recast as a 
prior distribution for an upcoming diagnosis to be made at tp + T. 

Completing the feedback loop from prognosis back to diagnosis in this manner has the potential to enhance 
the identifiability of damage in a model-based (inverse problem) SHM approach. A primary challenge of such 
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methods is that the inverse problem is typically ill-posed,!° e.g., a possibly infinite number of damage states 
can reproduce the sensor data. This problem is exacerbated with decreased signal-to-noise ratio resulting 
from sensor error, sparse measurements, small damage size, and limited sensor proximity to damage. In 
deterministic approaches to inverse problem SHM, a carefully-tuned regularization strategy must be employed 
during the solution process to ensure a unique solution exists.‘' Under a Bayesian regime, regularization 
is carried out more naturally through the specification of an informative prior distribution, favoring those 
solutions to the inverse problem that agree with any available a priori information, including expert opinion 
and results from previous analyses. While prior distributions for the location of damage initiation can be 
developed from offline hot spot identification, any reasonable prior distribution for damage size or severity 
must be generated online using prognostics. 

In a general SHM Bayesian inference problem, a parameter-dependent model response, f(Q), experimen- 
tal measurements, Y, and measurement errors, €, are related through the statistical model 


T=f(Q)+e. (1) 


Here, Q, ¢, and YT are random variables. The variable Q denotes the model parameters and has realizations q 
while Y has realizations in the form of observations, denoted by v. For diagnostics, the model response would 
be a directly measurable quantity like dynamic response or strains, and the parameters would characterize 
the damage in a given structure or component. For prognostics, the model response would be of the same 
form as the diagnostic output, and the parameters would be input to some predictive model (e.g., crack 
growth rate parameters). The goal of the Bayesian inverse problem in the context of SHM is to determine 
the posterior density, 7(q|vVons), that reflects the best or most likely distribution of parameter values based 
on the available observations. The multi-variate parameter distribution can be determined using Bayes’ 
Theorem, 


_ mviq)to(q) — st (v|q)t0(q) 
Naa). fuer elaymu ade 2) 


While the details will be omitted for brevity, obtaining a direct solution for the posterior in Equation 
2 is often intractable due to the evaluation of the normalizing constant in the denominator. However, 
sampling techniques exist!* that enable the approximation of the posterior given the likelihood distribution, 
m(ul|q), (determined through a comparison of model responses and observations) and the prior distribution, 
to(q), of the parameters. To avoid influencing the inverse problem with inaccurate prior knowledge, it 
is common practice to use a noninformative prior distribution; i.e., an improper uniform density over a 
bounded support.!?:!3 However, as mentioned above, there are significant benefits to an informative prior 
distribution. To reap these benefits, the prior distribution at a given time after a prognosis has been made, 
t=tp) +T, is 


T0(4, ty + T) = P(A, ty so T, D;,), (3) 


where P is the stochastic response of the predictive model extrapolated forward in time by 7 given a set of 
model parameters, #, and the set of damage diagnoses available at the time of prediction, D,,,. 


III. Numerical Example 


In this work, a numerical example is used to demonstrate the diagnosis-prognosis coupling. An alu- 
minum alloy 2024-T3 tensile specimen was modeled using finite elements,'* and the fracture mechanics code 
FRANC3D?° was used to insert a through-crack of length 0.0735 in into the model and simulate its growth 
under constant amplitude fatigue loading. The maximum applied stress was 6.25 ksi and the minimum 
applied stress was 0.625 ksi. The crack was grown using linear elastic fracture mechanics and a median ex- 
tension approach.!° An illustration of the specimen and associated boundary conditions is shown in Figure 
2. Additional information regarding the model can be found in the Appendix. 

Stress intensity factors (SIFs), denoted as K, were computed via finite element analysis of the cracked 
mesh and used to define the vector of crack extensions over the discretized crack front. The crack was 
extended at each step such that the median growth was equal to a user-defined value of 0.005 in. The crack 


growth rate, 44, was defined using the NASGRO (version 4) equation!® as 
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dt toe 


d 1- AK 

. c|( r)AK| a (4) 
with units of in/cycle. The parameters C, n, p, and q are empirical constants. The load ratio, R, is the ratio 
of the minimum SIF, Kinin, and maximum SIF, Kaz, for a given cycle and is time-invariant for constant 
amplitude loading. The Newman closure term is denoted by f, AK is the SIF range defined as Kingz —Kmin, 
AK, is the threshold SIF range below which no growth occurs, and Kc is the fracture toughness. The 
parameter values used in the present work can be found in the Appendix. The crack growth simulation was 
considered ground truth in place of experimental results and was used to provide simulated measurements 
and to quantify accuracy when investigating the proposed hypothesis. 
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Figure 2. Illustration of the numerical example. The left diagram shows the nominal dimensions of the specimen and 
the applied boundary conditions where o is the applied stress. The fixed and roller regions imply no displacement 
and displacement only in the y direction, respectively. The crack length, a, changes over time and is measured via a 
strain-based diagnostic technique.’” Simulated strain sensors are placed a distance h from the crack, which is located 
at the mid-plane of the specimen. The right side of the figure is a qualitative illustration of the strain field evolution as 
the crack grows from the first measurement to failure. The entire simulation consisted of 250 median extension steps. 


Throughout the specimen lifespan, a strain-based damage diagnosis technique!” based on Bayesian infer- 
ence was used to estimate the crack length with quantified uncertainty (i.e., g = a in Equation 2). Strains 
were extracted as shown in Figure 2 along two rows of simulated sensors, parallel to the direction of crack 
growth and at a distance, h = 0.93 in. For the purposes of this study, the simulated measurements were aug- 
mented with Gaussian white noise according to a signal to noise ratio, SNR = 20. The prognosis framework 
was a simplified version of an approach!® that utilized Bayesian inference and the same NASGRO (version 
4) crack growth rate equation (Equation 4) used in the FE-based simulation. In contrast, however, the crack 
growth rate equation was solved directly as a one-dimensional approximation of the three-dimensional crack 
growth model. To improve this approximation and reduce sources of error that might interfere with the study, 
the mean SIF's obtained from each crack step of the FE-based simulation were placed in a lookup table. A 
linear interpolator was then used when numerically solving the ordinary differential equation (Equation 4), 
which is a function of a and t, noting that AK is a function of a, R, and o. 

A set of diagnoses was required in order to obtain the first prognosis and initialize the diagnosis-prognosis 
feedback loop. Throughout the crack growth simulation, diagnoses were made at the fixed median extension 
intervals defined during the fracture simulation and corresponding to the set Tp = {t3,...,tp} where tp = ta7. 
Times t; and tz were discarded. For simplicity, each probabilistic diagnosis was dimensionally reduced to a 
single mean crack length in this study and placed in the set D;,. At time tp, D;, was used to predict the 
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distribution of crack lengths for the set of times Tp = {t2g,...,t39}. These distributions were passed to the 
diagnostic framework as prior knowledge for the subsequent diagnoses at all times t € Jp. The diagnosis 
results with and without the prognosis feedback were investigated to quantify the effectiveness of the method. 


IV. Results & Discussion 


Figure 3 shows the prognosis results generated for 7p U Jp based on diagnoses from 7p only. The 
prognosis output consists of a mean prediction (dashed black line) and the 95% credible and prediction 
intervals, all of which can be compared to the ground truth curve obtained from the numerical crack growth 
simulation. The diagnosis data (i.e., mean crack lengths) obtained for all t € Tp are plotted as white circles. 
At times t € Jp, pairs of diagnoses were made corresponding to those with and without feedback. These 
diagnoses retained their full distributions, but only the means are plotted in Figure 3(a) for clarity. Figure 
3(b) zooms in on Tp and represents the diagnosed crack length distributions as error bars corresponding to 
the bounds that contain ~ 95% of the distribution’s probability. 
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Figure 3. Prognosis and diagnosis results overlain on the ground truth crack growth curve (solid black line). Diagnoses 
up to time t, used to generate the prognosis are represented by white circles. The prognosis includes a mean prediction 
(dashed line) as well as 95% credible and 95% prediction intervals. Sub-figure (a) only shows the mean crack lengths 
generated by the diagnosis framework for clarity. Sub-figure (b) is a zoomed view of the times after t, at which 
feedback was used and includes the 95% probability bounds on the diagnosis distributions to demonstrate the effect of 
the feedback on uncertainty. 


From Figure 3(a), the prognosis is accurate in the short term, evidenced by the fact that the true crack 
lengths all lie within the credible interval for Tp. The prediction interval encompasses ~ 95% of the data, 
as expected. Additionally, the diagnoses with feedback from this prognosis are generally more accurate 
than those without as they tend to be focused more around the true crack length than the diagnoses with 
noninformative prior distributions. Examining Figure 3(b), the effect on the means becomes clearer. At all 
times, the prognosis-supplied prior distribution tends to pull the mean of the uninformed diagnosis toward 
the prognosis mean value. As a result, if the uninformed diagnosis mean lies between the prognosis mean 
and the true crack length for a given time, the resulting informed diagnosis mean is less accurate, albeit by 
a small amount since the prior and the likelihood distributions are similar in these cases. The interaction 
between the noninformative (i.e., the likelihood) and the informative prior distributions is illustrated more 


5 of 9 


American Institute of Aeronautics and Astronautics 


clearly when examining the full marginal probability density functions (PDF) of the diagnoses and prognosis 
at times t € Jp, as shown in Figure 4. For all values of r investigated, when the uninformed diagnosis mean 
lies outside the region bounded by the prognosis mean and the true value, accuracy is increased. 
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Figure 4. Probability density functions (PDF) for the first six diagnoses made at t = (a) 577,650, (b) 584,848, (c) 
591,784.0, (d) 598,410, (e) 604,684 and (f) 610,812 cycles. For each time, the uninformed diagnosis and feedback- 
informed diagnosis are shown along with the prognosis made at t,, which was used to generate prior distributions for 
the feedback case. 


Accuracy was only the first metric to be investigated; a reduction in diagnosis uncertainty was also 
desirable. Figure 3(b) and Figure 4 both show that the uncertainty is decreased in all cases where feedback 
is incorporated. This is indicated in the former by smaller error bars, meaning that the region of 95% 
probability is more concentrated than is the case for the uninformed diagnosis. In the latter figure, the 
reduction in uncertainty is represented by a decrease in width and increase in height of the crack length 
PDFs. All of these comparisons are quantified in Table 1 along with measures of accuracy. Here, the percent 
error between the true crack length, atrue(t), and the mean diagnosis for all t € Tp was calculated as 


sal ° 


where d(t) is the mean crack length obtained from the diagnosis at time t. Similarly, the reduction in the 
width of the 95% probability bounds was calculated as 


% Error = 


wp(t) — wni(t) (6) 
WNI (t) : 
where wp(t) and wyy(t) indicate the width of the 95% probability bounds for the feedback and uninformed 
cases, respectively. For the present example, the addition of prognosis feedback resulted in a decrease in 
diagnosis uncertainty. Accuracy was generally improved except for the special cases where the diagnosis 
was already in fairly good agreement with the prior distribution provided from the prognosis. For cases 
where a decrease in accuracy was observed, the decrease was negligible and could be an artifact of sample 
size. Despite the promising performance shown, the results are necessary but not sufficient to prove the 
diagnosis-prognosis feedback hypothesis. Care was taken to isolate the effect of the feedback, but more work 
will be required to ensure that the diagnosis improvement was not based on other factors. 
Furthermore, it is important to address the potential re-use of data, which would violate tenets of the 
underlying Bayesian analysis. It was assumed herein that, since the information from prognosis was only 


Change = 
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Table 1. Quantitative effect of the stochastic feedback loop. For the set 7p, the percent error between the mean crack 
length provided by the diagnoses and the true crack length is given for the noninformative prior distributions and 
informative prior distributions (i.e., from feedback), respectively. The better performance at each time is highlighted 
with bold font. The percent change in the spread of the diagnosis distribution induced by the inclusion of the prognosis 
feedback relative to the diagnoses without feedback is also presented. Here, a negative sign means that the uncertainty 
has decreased via feedback. 


Time relative Percent error in mean Percent error in mean Percent change in 95% 
to ty diagnosis (noninformative) | diagnosis (feedback) probability bound width 
ty4+1 16.70 8.60 -25.66 
ty+2 -1.95 -0.82 -24.97 
tors 8.81 4.64 -36.43 
ty+4 -1.47 -0.09 -28.37 
tr+s -3.92 -0.95 -29.87 
ty+6 0.39 1.24 -22.86 
tp+7 -5.52 -1.37 -20.91 
to18 -4.63 -2.36 -16.97 
ty+9 3.25 3.58 -29.68 
tp+10 -5.42 -1.61 -12.90 
tp4i1 -3.71 -1.36 -12.30 
a 2.69 3.50 -17.28 


reused in the form of prior knowledge, there was little difference between the automated feedback and any 
given expert opinion, which is a typical source of prior information. Additionally, the original data gathered 
from the numerical simulation were strains, and there is no re-use of these strain values except what is 
processed and propagated forward through the prognostic framework. The feedback itself, in this case, 
consisted of a crack length distribution and no stress or strain information. The authors argue that this is 
enough, especially when coupled with the demonstrated improvements in both accuracy and uncertainty, to 
warrant further investigation into the approach. 


V. Concluding Remarks 


Traditionally, the relationship between damage diagnosis and prognosis is unidirectional. An alternative 
approach was investigated herein where diagnosis and prognosis are instead related through a stochastic 
feedback loop. Prognosis in the structural health management or digital twin context involves simulating 
damage evolution from the current state to failure. If probabilistic in nature, model output in the form 
of damage state distributions are available at the times between the current state and failure. These times 
correspond to future diagnoses, and, as such, these predictive damage state distributions can be used as prior 
distributions in a Bayesian-based diagnostic framework. To test the feasibility of this concept, a numerical 
fatigue crack growth simulation was conducted. The diagnosis-prognosis feedback loop was implemented, 
and diagnosis results with and without the prognosis-generated feedback were compared. In general and for 
all times tested, the feedback loop resulted in a sizable reduction in uncertainty as well as either an increase 
in accuracy or a negligible reduction in accuracy (i.e., no significant change). 

The results obtained in the study are promising and support future research into the concept of a 
diagnosis-prognosis feedback loop. Many avenues remain open to investigation, and future work should 
include a deeper look into whether or not the potential re-use of data violates any aspect of the underlying 
Bayesian analyses. Care was taken in this study to limit noise levels and ensure a quality prognosis, resulting 
in an ideal case. The next step in this research will be to begin searching for degenerate cases where the 
feedback approach not only fails to improve the results but also degrades them (i.e., significantly decreases 
accuracy or increases uncertainty). While not necessarily classified as a degenerate case, there may be a 
time horizon beyond which prognosis results are no longer beneficial depending on the uncertainty in the 
prediction; this problem could potentially be alleviated by continuously updating the prognosis as often as 
possible. More importantly, as with any feedback loop, divergence is an issue that should not be taken 
lightly. Cases which cause instability or divergence should be identified if they exist. Finally, other forms of 
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the diagnosis-prognosis feedback loop should be investigated, such as those that could be applied directly to 
the measurement process. Defining the limits of the proposed feedback approach as well as exploring new 
ways in which it can be applied could yield significant benefits for both structural health management and 
digital twin. 


Appendix 


The results presented herein were dependent on a numerical crack growth example that served as the 
ground truth for gathering simulated data and evaluating accuracy of the proposed approach. The initial 
finite element model mesh used for this simulation comprised 3,699 Tet10 elements and 7,024 total nodes. The 
analysis was completed using the implicit version of the parallel finite element code ScIFEN.'* Remeshing 
was conducted within FRANC3D?° upon crack insertion and at every growth step. The mesh size was then 
dependent on the length of the crack and the crack front template radius, which is a special ring of elements 
placed around the crack front. In this study, the template radius was 0.002 in. The elastic modulus and 
Poisson’s ratio used for aluminum alloy 2024-T3 were 10,600 ksi and 0.33, respectively. 

The NASGRO parameters used to model the crack growth rate in the numerical simulation are presented 
in Table 2. The NASGRO equation (Equation 4) can be fully defined with these parameters, and is a function 
of only crack length, a, and time, t, measured in fatigue cycles. The dependence on a is introduced via the 
stress intensity factor (SIF) range, AK where AK = Kaz — Kmin. Kmaz and Kmin are the maximum and 
minimum SIF's at each cycle, respectively, and depend on the crack length, a, and the applied load. Some 
of the parameters are functions of more fundamental material properties; in that case, the equations used 
to calculate them are provided. The load ratio, R, is the ratio of the minimum and maximum SIF for a 


Table 2. Parameters used in the NASGRO equation (Equation 4). 


Parameter Value Equation 

C 8.0x 10-9 - 

n 3.20 - 

p 0.25 - 

q 1.00 - 

R 0.10 pains Geter 

AKin ‘ AK$ [2B] FRC 701 — Ay) -RICK, 
AK? : AK, [2]? 

ag 1.5x10-3in - 

ane 1.21 : 
AK, 1.22 ksivin  - 

f - max(R, Ag + AiR + A2R? + A3R?) 
Ao : (0.825 — 0.340 + 0.0502) x [cos (% Smaz)]/* 
Aj z (0.415 — 0.071@) ma= 
As ‘ be Ag eAi— As 
As Z 9Ag Ay A 

Smaz 0.30 . 
o0 
a 2.00 - 
Ke - Kic(1 + Bye Axt/to)” 
Kic 30.0 ksiVin = - 
Ax 1.00 : 
By 1.50 2 

to : 2.5(Krc¢/oy)? 
oy 53.0 ksi - 

t 0.08 in Z 
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given cycle (R is constant for the simulations presented herein), f is the Newman closure term, AKy,, is the 
threshold SIF range below which no growth occurs, AK is the observed threshhold for a high R value, Ko 
is the fracture toughness, Kyc is the plane strain fracture toughness, a is a plane stress/strain constraint 
factor, Smax is the ratio of maximum applied stress to flow stress, t is the approximate thickness of the 
specimen, oy is the material yield stress, and ap is a small crack parameter. The remaining parameters are 
empirical constants. 
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